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This paper examines the relation between the spin- wave instabilities of collinear magnetic phases 
and the resulting non-collinear phases for a geometrically- frustrated triangular-lattice antiferromag- 
net in the high spin limit. Using a combination of phenomenological and Monte- Carlo techniques, 
we demonstrate that the instability wave- vector with the strongest intensity in the collinear phase 
determines the wave- vector of a cycloid or the dominant elastic peak of a more complex non-collinear 
phase. Our results are related to the observed multi-ferroic phase of Al-doped CuFeCb. 
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It is well-known that the transition between different 
magnetic ground states may be signaled by the soften- 
ing of a spin- wave (SW) mode. In the simplest case of 
a conventional square-lattice antiferromagnet, the soft- 
ening of a SW mode at wave- vectors (-zr, 0) and (0, 7r) 
signals the spin flop and canting of the magnetic mo- 
ments at a critical field. In the manganites the SW 
instabilities of the ferromagnetic state have been used 
to construct the phase diagram for the antiferromagnetic 
(AF) phases that appear with Sr doping. The soften- 
ing of a SW excitation at (7r, 0) signals the instability 
of the Neel state and the canting of the spins in a spin- 
1/2 union-jack lattice Q. But the relation between the 
SW instabilities of a collinear phase and a resulting non- 
collinear phase is less clear when multiple SW instabilities 
occur simultaneously or when the non-collinear phase has 
a complex magnetic structure with several elastic peaks. 
This paper explores the relation between the SW instabil- 
ities of the collinear 4 and 8-sublattice (SL) phases of a 
geometrically-frustrated triangular-lattice antiferromag- 
net (TLA) and the non-collinear phases that appear with 
decreasing anisotropy D. When multiple SW instabilities 
of the collinear phase occur at once, the instability wave- 
vector with the largest intensity determines the dominant 
ordering wave- vector of the resulting non-collinear phase. 
One of the predicted non-collinear phases may be related 
to the multi-ferroic phase that appears in CuFe02 with 
Al doping 

Frustrated TLAs with AF nearest-neighbor exchange 
Ji < exhibit a remarkable number of competing ground 
states [4] . With interactions Ji up to third nearest neigh- 
bors (denoted in Fig. [1]) and assuming Ising spins along 
the z direction, Takagi and Mekata [5] obtained a phase 
diagram with ferromagnetic (FM), 2-SL, 3-SL, 4-SL, and 
8-SL phases. A portion of that phase diagram is sketched 
in Fig. [TJ For the geometrically-frustrated TLA CuFe02 
in fields below 7 T, the ground state is the 4-SL phase 
@, Q and the black dot in Fig. [T] denotes the esti- 
mated ratio of exchange parameters J2/I Ji| ~ —0.44 and 
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FIG. 1: (Color online) A portion of the TLA phase diagram 
for large D and J\ < 0. The dashed (red) curve divides the 
4-SL phase into 4-SL I and 4-SL II regions. For the 4 and 8- 
SL phases, red circles denote up spins and blue circles denote 
down spins; the solid (gray) line defines the unit cell. The 
dashed (blue) line obeys the relation J3 = I.3J2, the black 
circle is the estimated location of the exchange parameters for 
CuFeCb, and the white circle lies on the boundary between 
the 4-SL I and 4-SL II regions. 



J 3 /|Ji|»-0.57 1,0. 

As demonstrated by the small SW gap of about 0.9 
meV 0, [H on either side of the ordering wave-vector 
q = 7rx, the spin fluctuations of CuFe02 are much softer 
than would be expected for Ising spins. With Heisen- 
berg spins, the collinear magnetic phases of a TLA be- 
come locally unstable below a critical anisotropy D c that 
depends on the exchange parameters Ji. The observed 
softening of the SW modes in CuFei-^ALcC^ with Al 
doping [3[ can be reproduced by lowering D towards D c 
|l0[ in the 4-SL I region of Fig. [TJ For Al concen- 
trations above about x c « 0.016, the magnetic ground 
state of CuFei_ x Al a; 02 becomes non-collinear and dis- 
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FIG. 2: (Color online) Energy as a function of D/\J\\ for 
the 4-SL (black dashed), CNC (blue diamonds), and cycloid I 
(red circles) phases with J3/J2 = 1.3 and J1/J2 = 2.28. The 
bottom diagram shows the the cycloid for arbitrary q. The 
top diagram shows the CNC phase with up spins in red and 
down spins in blue. 



plays multi-ferroic properties pjj, LL2|, [13| . 

We have determined the ground-state magnetic phases 
of the TLA using a combination of Monte-Carlo (MC) 
simulations and phenomenological techniques. The TLA 
Hamiltonian is 



H 



S 7 ; • S 7 ' 



(i) 



where Jij includes first, second, and third- neighbor inter- 
actions (shown in Fig. [1]). The nearest-neighbor distance 
has been set to 1. The SW frequencies of the collinear 
phases are obtained by performing a Holstein-Primakoff 
1/5 expansion about the classical limit. Above the crit- 
ical anisotropy D C1 a collinear phase is locally stable if 
the SW frequencies o;(k) are positive and real for every 
momentum k. 

MC simulations were used to find the non-collinear 
magnetic phases of the TLA. The simulations were 
started at a high-enough temperature to rule out 
metastable states. To mimic the process of thermal an- 
nealing, the system was slowly cooled to a final tempera- 
ture (in units of | J\\S 2 ) ranging from 4x 10 -3 to 1 x 10 -4 . 
Lowering the final temperature further did not signifi- 
cantly change the resulting non-collinear phase. Using 
lattices of varying sizes with periodic boundary condi- 
tions, we found that there was no substantial change for 
lattices greater than 16 x 16. 

In Fig. [TJ the 4-SL phase is separated into regions I 
and II by the curve J\j J 2 — 2 = J 2 / J3. In region 4-SL II, 
the instability wave- vectors are given by k = (7r±7r/3)x, 
independent of the exchange parameters; in region 4-SL 
I, the instability wave-vectors depend on the exchange 



parameters [14j. With the exchange parameters corre- 
sponding to the black circle in region 4-SL I, we deter- 
mined the stable magnetic phases as a function of D. 
As shown in Fig. O the 4-SL phase is stable down to 
D/\J\\ w 0.27, below which MC simulations obtain the 
complex non-collinear (CNC) phase shown on the top of 
Fig. O Since the same up or down spin frequently occurs 
at sites R = rax + n\/3y and R/ = R + x/2 + \/3y/2 or 
R" = R — x/2 + v^y/2, the CNC phase retains some of 
the FM correlations present in the 4-SL phase. Transla- 
tionally invariant in the y direction, the CNC phase has 
a period in the x direction between 2 and 3 lattice con- 
stants. Becauase the MC simulations were performed on 
a finite lattice, the energy of the CNC phase is overesti- 
mated and we cannot say whether this phase is commen- 
surate or incommensurate. 

Below a second threshold value of D/\J\\ w 0.08, a 
cycloid like the one sketched in the bottom panel of Fig. [2] 
[15[ has a lower energy than the CNC phase. As discussed 
below, the wave-vector of the cycloid is independent of 
D. If the CNC phase were neglected, then the cycloid 
would achieve a lower energy than the 4-SL phase below 
D/\J\\ ^0.2, still above the critical value D C /\J\ \ ^0.15 
for the local stability of the 4-SL phase. 

To gain a better understanding of the phases stabilized 
within the TLA, we have evaluated the magnetic phases 
along the line with J3/J2 = 1.3 drawn through the black 
dot in Fig. [U Five stable phases are presented as a 
function of \Ji\/D and \J 2 \/D in Fig. El The 4-SL phase 
is stable along a strip through the diagonal of this plot. 
Although not indicated by this figure, the 4-SL region 
disappears above \Ji\/D w 40. Close to the origin or 
for large D, a collinear 8-SL region is indicated in Fig. 
[U Two cycloids are also obtained: in the upper left, 
cycloid II with wave- vector 47rx/3; in the lower right, 
cycloid I with the variable wave- vector indicated in the 
figure [Hj]. Finally, a CNC phase appears just below the 
4-SL phase and disappears above \J\\/D w 20. Regions 
of local stability for the collinear phases are indicated in 
Fig. [3] by the dashed black lines [14[. The results in Fig. 
[2] can be obtained from Fig. [3] by drawing a line from the 
origin with slope 2.28 (gray) (so that J 2 /\J\\ = —0.44), 
which passes from the 4-SL phase through the CNC phase 
into cycloid I. 

The classical energies of each of these phases can be 
written as E/S 2 = A 1 J 1 + A 2 J 2 + A 3 J 3 - A D D. The 
coefficients for each phase are given in Table HI Only a 
non-collinear phase with 0.5 < Ajj < 1 can intercede 
between a collinear phase with Ad = 1 and a cycloid 
with A D = 0.5. For the CNC phase with A D w 0.71, the 
error bars indicate the range of parameters obtained from 
MC simulations near \Ji\/D = 5.7 and \J 2 \/D = 2.5. 
This phase is characterized by rather weak next-neighbor 
correlations with small \A 2 \. The CNC phase space in 
Fig. [3] is obtained by using the results of Table U to 
evaluate the energies of the MC spin configurations as 
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gion can be related by a symmetry operation to one of 
these three. We find that the instability at ki always 
has a larger intensity than the "twins" at k2 or k3 or 
than any of the other wave- vectors related by symmetry. 
Correspondingly, cycloid I along any diagonal in Fig. [3] 
has the same wave-vector q as the dominant instability 
of the 4-SL phase. 

Similar conclusions are reached for the 8-SL phase, 
which switches to cyloid I along any diagonal in Fig. 
02 Although the SW instability of the 8-SL phase oc- 
curs simultaneously at several wave- vectors, the domi- 
nant wave- vector instability of the 8-SL phase coincides 
with the wave- vector q of cycloid I along any diagonal in 
Fig. El 

However, the CNC phase that intercedes between the 



FIG. 3: (Color online) Phase diagram for the TLA as a func- 
tion of \Ji\/D and \Ji\/D with J3 = I.3J2 containing five 
regions: 4-SL (blue), 8-SL (green), CNC (violet), cycloid I 
(variable green-orange), and cycloid II (maroon). The dashed 
(white) line separates regions 4-SL I and 4-SL II. The dotted 
(black) curves denote the metastable boundaries for the 4-SL 
and 8-SL regions. Cycloid I has wave-vectors q that range 
from 0.684tt to 0.923tt in intervals of 0.016tt. 



functions of \J\\/D and \J2\/D. Hence, the CNC region 
may be underestimated. 

The ordering wave- vector q = gx of cycloid I is eval- 
uated by minimizing E with respect to q. So q depends 
only on the ratios J2/J1 and J3/J2, as indicated by the 
diagonal lines in Fig. 03 Cycloid II with q = 47r/3 cor- 
responds to the 120° Neel state found in a classical TLA 
with D = and nearest and next-nearest neighbor inter- 
act ions|l6| when | J 2 / J\ | < 1/8. A slightly distorted Neel 
state [15| is stable for nonzero D over a range of exchange 
parameters with | J3/J2I > 1/2, so that the diagonal line 
in Fig. [T] passes through the 4-SL and 8-SL phases. MC 
simulations were used to confirm the stability of cycloids 
I and II in Fig. H 

With decreasing D or moving away from the origin of 
Fig. [3] along a diagonal, the 4-SL phase becomes unstable 
either to cycloid II or to the CNC phase. The white line 
bisecting region the 4-SL strip in Fig. [3] corresponds to 
the white point in Fig. [I] at the border between the 4-SL 
I and 4-SL II regions with J2/IJ1I = —0.36. In region 
4-SL II or above the white diagonal line, the 4-SL phase 
has instabilities at the wave- vectors (7r±7r/3)x. The SW 
intensity at the larger of these two wave- vectors always 
dominates and the 4-SL phase evolves into cycloid II with 
wavector 47rx/3. 

In region 4-SL I, the 4-SL phase has three unique SW 
instabilities: one at wave-vector ki along the x axis, 
another at k2 rotated by 7r/3, and a third at k3 ro- 
tated by — 7r/3. All three have the same magnitude with 
it/ 2 < ki < 7r. Other SW instabilities in the 4-SL I re- 
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FIG. 4: (Color online) (a-b) SW frequency and intensity ver- 
sus wave- vector k x for the 4-SL phase with \J\\/D — 5.5, 
J3 = I.3J2, where \J2\I~D varies from 1.30 to 2.67. (c) 
Fourier transform for the S x , S y , and S z components of the 
CNC phase with the same exchange parameters as above and 
|J 2 |/L> = 2.5. 



TABLE I: Energy Coefficients for Collinear, Cycloid, and CNC Phases 



Phase Ai A 2 A 3 A D 

TSL 1 1 1 1 

8-SL 111 

Cycloid I -(cos(g) + 2cos(<?/2)) -(l + 2 cos(3<?/2)) - ( cos(2<?) + 2 cos(g)) 1/2 

Cycloid II 3/2 -3 3/2 1/2 

CNC 0.595±0.001 -0.097±0.001 1.161±0.001 0.712±0.001 



4-SL and cycloid I phases is characterized by several elas- 
tic peaks shown in Fig. [4|(c). Within the precision of our 
MC simulations, the dominant wave- vector k x ~ 0.877T of 
the CNC phase coincides with the dominant instability 
wave- vector of the 4-SL phase that preceeds it. 

To demonstrate how the magnetic ground state evolves 
from cycloid II into the 4-SL phase and then into cycloid 
I, we plot in Fig. Ufa) and (b) the SW frequencies and 
intensities versus wave- vector for the 4-SL phase with 
J\/D = —5.5, J3/J2 = 1.3, and \J2\lD varying from 
1.30 to 2.67. As I J2\/D approaches the lower limit for the 
stability of the 4-SL phase, the SW intensity dominates 
at the cycloid II wave- vector q = 47rx/3. At the upper 
limit, the SW intensity dominates at the cylcoid I wave- 
vector q « 0.837TX. Hence, the wave- vector of the SW 
instabilities for the collinear phases correspond to the 
ordering wave- vectors of the non-collinear phases. 

Bear in mind that the transition between magnetic 
ground states is not always signaled by the softening of a 
SW mode. In fact, the transition between the 4-SL and 
8-SL phases in Fig. [1] occurs even at D = 00, when the 
spins are Ising variables and the SW gaps for both phases 
diverge. 

The CNC phase may be related to the multi-ferroic 
phase observed in Al-doped CuFe02 [H, which was re- 
cently investigated by Nakajima et al. [l7j |. Based 
on neutron-scattering measurements, those authors con- 
cluded that the ground state is a modified cycloid with 
the same spin on sites R and R/ (see above). This 
phase has peaks at wave-vectors on either side of 7rx, 
in agreement with the neutron measurements. However, 
a modified cycloid cannot be stabilized by a Hamiltonian 
with the form of Eq. (pQ) , regardless of the exchange and 
anisotropy parameters. With an additional phase slip S 
for the spins at sites R', a pure cycloid with S = and a 
single elastic peak always has lower energy than the phase 
proposed in Ref.[l3] with 5 = —q/2. This conclusion has 
been verified by MC simulations. 

Like the non-collinear phase proposed earlier [l7|, the 
CNC phase also contains FM correlations between sites 
R and R' or R". So the CNC phase also has elastic 
peaks on either side of 7rx at k x ~ 0.877T and 1.137T, as 
shown in Fig. Htc). Because the FM correlations are not 
perfect and vary along the x direction, the CNC phase 
contains several other elastic peaks that may allow it to 
be experimentally distinguished from the phase proposed 



in Ref.0. 

To summarize, we have shown that the dominant wave- 
vector of a non-collinear phase in a frustrated TLA cor- 
responds to the dominant instability wave-vector of a 
collinear phase as the anisotropy is lowered and spin fluc- 
tuations become softer. The CNC phase sketched in Fig. 
[2] is a more reasonable candidate for the multi-ferroic 
phase observed in Al-doped CuFe02 than the one previ- 
ously proposed. 
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